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ABSTRACT 

The possible role of magneto-rotational instability (MRI) and its driven MHD turbulence in 
the solar interior is studied on the basis of the linear and nonlinear theories coupling with 
physical parameters, providing solar rotation profile inverted from the helioseismic observa- 
tion and a standard model for the internal structure of the sun. We find that the MRI venue 
is confined to the higher-latitude tachocline and lower-latitude near-surface shear layer. It is 
especially interesting that the MRI-active region around the tachocline closely overlaps with 
the area indicating a steep entropy rise which is required from the thermal wind balance in 
the sun. This suggests that the MRI-driven turbulence plays a crucial role in maintaining the 
thermal wind balance in the sun via the exceptional turbulent heating and equatorward angu- 
lar momentum transports. The warm pole existing around the tachocline might be a natural 
outcome of the turbulent activities energized by the MRI. 
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1 INTRODUCTION 

Two important unsolved issues in solar physics are the mechanisms 
of magnetic dynamo and angular momentum transport in the solar 
interior. These two processes essentially affect one another non- 
linearly, which makes the physical processes more difficult to re- 
solve even using huge computer facilities. Despite the development 
of numerical and observational methods, we have not yet drawn a 
complete physical picture of these processes. 

Solar physicists believe that two types of dynamo mechanisms 
operate in the solar interior. The flux transport-type dynamo sus- 
tains the global, coherent, magnetic field by the stretching £2-effect 
around the tachocline (Parker 1993; Dikpati & Charbonneu 1999; 
Dikpati & Oilman 2006). The convective dynamo, in which the 
convective motion stores free energy, plays a role in generating rel- 
atively weak, small scale magnetic fields (Parker 1955; Weiss 1994; 
Miesch et al. 2000; Brun & Toomre 2002; Brun et al. 2004). 

The solar rotation profile revealed by helioseismic observa- 
tion tell us that our knowledge of angular momentum transport in 
the solar interior is still insufficient (Kosovichev 1997; Schou et al. 
1998; Thompson et al. 2003; Howe 2009). We have no consensus 
on how the differential rotation is physically developed and main- 
tained in the convective zone (e.g., Kitchatinov & Rudiger 1995; 
Durney 1999; Rempel 2005; Miesch et al. 2006; Balbus 2009; Bal- 
bus et al. 2009). 

There is a crucial MHD process related to these outstanding is- 
sues, which is "Magneto-Rotational Instability (MRI)" (Balbus & 
Hawley 1991, 1994, Menou et al. 2004; Parfrey & Menou 2007). It 
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is essentially of local natures and destabilizes the differentially ro- 
tating magnetized system with a radially negative shear rate, that is 
g = d In f2/d In cr < 0, where fl is the angular velocity, and m is the 
cylindrical radius. The typical spatial scale of unstable MRI mode 
inversely cascades with its linear growth, and it finally operates the 
MHD turbulence. The MRI-driven turbulence is broadly investi- 
gated as a promising candidate for angular momentum transporter 
in the astrophysical accretion disks (Balbus & Hawley 1998). 

The stability of the solar interior to the MRI is studied at first 
by Balbus & Hawley (1994) to understand why the solar radiative 
core rotates uniformly. Menou et al. (2004) construct a linear the- 
ory of the MRI in the general form which is more suitable for the 
stellar interior (including heat, magnetic, and viscous diffusivities). 
Recently, Parfrey & Menou (2007) have pointed out the possibility 
that the MRI-driven turbulence prevents the global coherent dy- 
namo activity in the MRI active higher-latitude tachocline. 

Here we investigate, on the basis of the linear theory, where 
the MRI should be active in the convective zone and tachocline by 
taking account of the observed rotation profile inverted from the 
helioseismic data and the standard solar structure. Then we discuss 
a role of the MRI-driven turbulence for maintaining thermal wind 
balance which would be achieved in the solar interior. 



2 MRI VENUE IN THE SOLAR INTERIOR 

We indicate the promising venue for the MRI in the solar inte- 
rior. Considering axisymmetric WKB plane wave perturbation with 
e.Tcrt[i(k m m + k,z - cot)] in the Boussinesq approximation, the sta- 
bility of the system to the MRI is governed by the following local 
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Figure 1. The maximum growth rate of the MRI y max which is derived from 
the dispersion equation (T} is demonstrated by the color map. Normalization 
is the local angular velocity £2. The region with the positive growth rate is 
the promising MRI venue. The rotation profile is overplotted by the solid 
contour line. It increase by 15 nHz from 330 nHz to 480 nHz. Note that the 
region with zero growth rate of the MRI is filled with white. 



dispersion equation, in the cylindrical coordinate (m, (f>, z), 



k 2 

K ~ 4 - 2 

— 0} m 0J K + 6J„„ 



— (Dp)D\nPp- y 
TP 



where 



: 0)„0J V - (k • V A ) 2 



■ 4£2 2 (k • \ A fw K = , 



(1) 



-oj] - (k ■ v A ) 2 



: w + i\k 2 (x = k, rj, v) , D = 



k w d d 
k, dz dm 



(see Menou et al. 2004 for details). Here k 2 = k 2 ^ + kl is the 
wavenumber, v A = B/(4^p) I/2 is the Alfven velocity, B is the mag- 
netic field, and y is the adiabatic index of the gas. The heat, mag- 
netic and viscous diffusivities are represented by k, tj and v respec- 
tively. The others have their usual meanings. 

Note that the dispersion equation {TJ is obtained under the 
condition with stratified hydrostatic background, not limited only 
to the convectively stable background. The convectively unstable 
mode is also included in this dispersion equation. The instability in 
the form of thermal convection is treated on a footing equal to that 
of rotational instability in this paper as was done by Balbus (2009). 
The strength of magnetic fields in the unperturbed state is assumed 
to be weak enough that it does not affect the hydrostatic balance. 

We can search the most rapidly growing MRI mode at an ar- 
bitrary meridian point (r, 6) in the solar interior by numerically 
solving the linear dispersion equation l[T) coupling with the ob- 
served rotation profile inverted from the helioseismic data of SOHO 
Michelson Doppler Imager (MDI) and the standard model for the 
solar internal structure (Christensen-Dalsgaard et al. 1996). When 
searching for the fastest growing modes, we limit their wavelength 
to being sufficiently smaller than the typical width of the tachocline. 

The promising venue for the MRI is demonstrated by the color 
map in Figure 1. The color scale represents the maximum growth 
rate of the MRI at each local meridian point. The vertical and hor- 



izontal axes denote polar and equatorial radii normalized by the 
solar radius. The rotation profile adopted in our analysis is over- 
plotted on the MRI map with solid contour which increases by 15 
nHz from 330 to 480 nHz. The angular velocity within 6 < 10° 
is extrapolated by the cubic spline method from the lower latitude 
helioseismic data, where 6 is co-latitude. 

The local structure and the strength of magnetic fields are as- 
sumed to be vertical and weak, that is fi- = 1 G here. The MRI 
venue is not drastically changed due to the field strength and struc- 
ture as long as we give a weak, vertical component of the magnetic 
field. The destabilizing effect arising from the unstable internal 
gravity wave is eliminated so as to better focus on the destabilizing 
effect of the MRI modes on the system. In contrast, the stabilizing 
effect due to the density stratification on the MRI is consistently 
introduced in our analysis. 

It is found that, in the solar interior, the MRI-active region 
is confined to the higher-latitude tachocline and the lower-latitude 
near-surface shear layer. The stabilizing effect due to the density 
stratification on unstable MRI modes is reduced significantly by 
double diffusive effect in the convectively stable tachocline region 
(Menou et al. 2004; Masada et al. 2007). The MRI venue thus cor- 
responds trivially to the differentially rotating zone with negative 
radial shear = d In fl/d In nr < here. In addition, the wavevector 
of the most unstable MRI mode tends to be parallel to the grav- 
ity vector in the convectively stable tachocline as is suggested by 
Balbus & Hawley (1994). This is because such modes can grow in 
decoupling with the stabilization effect of the density stratification. 
When focusing naively on the MRI mode with neglecting the effect 
of the convectively unstable mode, the MHD turbulence powered 
by the MRI would be being operated in these MRI-active regions. 

It would be natural from a physical perspective to discuss the 
impact of convective motion on the MRI. Since the convective 
mode (unstable mode of internal gravity wave) has a growth rate 
which is at least an order of magnitude larger than that of the MRI, 
the profile of large-scale differential rotation drawn using a helio- 
seismic time-averaging method should be modified by the vigorous 
convective motion and not be maintained during the typical evolu- 
tion time of the MRI. The stationary background assumed for de- 
riving the dispersion relation is thus not justified to be applied to the 
MRI in the sun's convective zone. When considering the vigorous 
convection, the growth of the MRI mode should be inhibited in the 
near-surface shear layer where is one of promising MRI venues. 

In contrast, the tachocline which is a thin transition region be- 
tween radiative and convection zones generally has a convectively 
stable nature with sub-adiabatic entropy gradient. The large-scale 
differential rotation should be thus maintained without being af- 
fected by the turbulent convection during the evolutionary phase 
of the MRI. The dispersion relation adopting the stationary back- 
ground as unperturbed state is applicable to the MRI mode (Par- 
frey & Menou 2007). In addition, as is denoted in this section, the 
stabilization effect due to the sub-adiabatic entropy gradient is re- 
duced by diffusive relaxation and does not have a crucial role in 
suppressing the growth of the MRI. We can thus expect that, in the 
tachocline region, the local scale MRI mode can evolve to the larger 
scale with leaving its character unchanged and finally operates the 
MHD turbulence with a typical size comparable to the radial width 
of the tachocline. The nonlinear interaction between the convective 
overshoot and the MRI-driven turbulence is discussed later in § 3. 

In the following, we focus on the role of the MRI-driven tur- 
bulence around the tachocline region. The impact of driven MHD 
turbulence on the thermal wind balance in the solar interior is a 
special interest of this paper. 
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Figure 2. (a) Specific entropy gradient in 6 direction obtained from the thermal balance equation (2) coupling with the observed rotation profile. The nor- 
malization is chosen as -10 Cp (S = S/Cp) here, (b) Temperature deviation AT from the equatorial plane calculated from the thermodynamic relation 
ST — TCp 1 J SS . The vertical and horizontal axes have same meanings as fig. 1 . The isorotation contour is overplotted just for the reference on both plots. 



3 POTENTIAL RELATION BETWEEN THERMAL WIND 
BALANCE AND MRI-POWERED TURBULENT 
ACTIVITY 

3.1 Entropy Rise Resulting from Thermal Wind Balance 

The effect of the MRI-driven turbulence on the global coherent dy- 
namo has already been argued by Parfrey & Menou (2007). Ac- 
cording to their logic, the coherent stretching process of the mag- 
netic field by the SI effect is prevented by the MRI-driven turbu- 
lence in the higher-latitude tachocline. This implies that the flux 
transport type dynamo efficiently works only in the lower-latitude 
tachocline. They insist that this prevention effect of the MRI-driven 
turbulence on the dynamo is why the active region is restricted 
within ~ ±45° around the equator on the solar surface. This is an 
interesting point of view for qualitatively explaining the formation 
mechanism of the activity belt. 

From a different physical viewpoint, we discuss here the role 
of the MRI-driven turbulence in maintaining the thermal wind bal- 
ance in the solar interior. The thermal wind equation with a baro- 
clinic term is 



dSl 1 

T 

dz 
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(Pedlosky 1987; Tassoul 2000) where 5 is the specific entropy, g 
is the gravity in hydrostatic equilibrium, and C P is the constant 
pressure specific heat. This follows from the time-steady vorticity 
equation in conditions with small Rossby numbers. Equation {2} in- 
dicates that the latitudinal entropy variation drives the differential 
rotation and yields the deviation from Taylor-Proudman state with 
dCl/dz = 0. 

When the observed rotation profile is given, we can directly 
calculate the left hand side of equation (O and derive the entropy 
gradient in 6 direction. Moreover, we can retread it to the tem- 
perature deviation by using the thermodynamic relation dT/88 ^ 
TCp'dS jdB because of the quasi-incompressible plasma (dP/89 ^ 
0), where T is the reference temperature varying only in the radial 
direction. This yields information about the latitudinal temperature 
inhomogeneity in the solar interior (e.g., Miesch 2005). 

Figure 2a illustrates the normalized latitudinal entropy gradi- 



ent -C p ' 10 6 35/30 consistent with the thermal wind balance. Fig- 
ure 2b demonstrates the latitudinal temperature deviation AT (= 
TCp 'EAS ). We choose the equatorial plane as a fiducial location, 
so AT(8 = 90°) = K in Figure 2b. Note that the temperature ob- 
tained from the standard solar structure (Christensen-Dalsgaard et 
al. 1996) is adopted as the reference temperature f at each radii. 
The isorotation contour is overplotted on both maps. 

The specific entropy is found to rise steeply toward the pole 
between 8 20° and 60° at around the tachocline depth. This sug- 
gests that an exceptional heat production is maintained in this area. 
In addition the temperature deviation AT increases with increasing 
latitude and reaches up to AT ~ 10 K at the polar region around the 
tachocline. We interpret physically this steep entropic rise and the 
temperature enhancement around the higher-latitude tachocline in 
relation to turbulent activities powered by the MRI in the following. 



3.2 MRI-Powered Turbulent Activity 

The higher-latitude tachocline where the entropy steeply rises is 
reminiscent of the MRI-active region denoted in Figure 1 although 
the governing equations providing each map are completely differ- 
ent. Then we attempt to reveal a potential relationship between the 
steep entropy rise and the MRI-powered activities around the solar 
tachocline. 

A nonlinear outcome of the MRI activity in astrophysical sys- 
tems is the enhancement of turbulent viscous heating. We can ex- 
pect that the viscous heating maintained by the MRI-driven tur- 
bulence plays an essential role in inducing the exceptional heat 
production around the higher latitude tachocline which is required 
from the thermal wind balance. 

In addition the angular momentum tends to be transported 
equatorward by the MRI-driven turbulence in the stably strati- 
fied tachocline (Balbus & Hawley 1994). Then the gaseous matter 
would advect toward the pole with absorbing the heat generated by 
the MRI-driven turbulence. The warm pole is thus expected to be 
formed at around the tachocline depth as a natural consequence of 
the MRI-powered turbulent activities. 

For simplicity, we assume that the differential rotation in the 
tachocline is maintained sufficiently long to balance the angular 
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momentum transport entirely in the solar interior. Then we can 
evaluate the turbulent heating rate provided by the MRI-driven tur- 
bulence in the tachocline region. 

The Maxwell stress resulted from the MRI-driven tur- 
bulence would tap into the fraction of the free energy stored in the 
differential rotation e diff and is given by 



M md> = — 



B m B,i, 



An 



= few ■ 



(3) 



where / is an arbitrary parameter ranged from 0.01 to 0. 1 (Masada 
et al. in prep). The free energy stored in the differential rotation e d iff 
would be written by 



ediff = Po [Q. z m (m)m z - Sl^Tj^dmdfidz / dmd<pdz 
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where 



frigid = [(wo + Am/2)/m Q ] cl . 

Here fldiff(ra-) is the angular velocity at differentially rotating state 
with maximum free energy and Q.n g n shows the constant angular 
velocity at rigidly rotating state with minimum free energy. Note 
that tjq is the fiducial radius, flo is the angular velocity at the fidu- 
cial radius, Am is the radial width of shear layer, and p is the av- 
eraged density around ot . 

The viscous heating rate resulted from the MRI-driven turbu- 
lence is thus, at around the fiducial radius C7 , 
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\q\ 2 Q.lzj Azj . 



where f = M^/\q\Clo is the viscous energy deposition rate. Note 
that the turbulent viscosity is then given by v tllrb - £/p. 

Figure 3 shows the viscous heating rate e which would be 
sustained by the MRI-driven turbulence around the tachocline as 
a function of co-latitude 6 in the cases / = 0.01, 0.05, and 0.1. 
The fiducial parameters Clo, q, and Act are systematically de- 
termined from the observed rotation profile by averaging over their 
radial variations. Thus, they vary only with the latitude. The av- 
eraged density around the tachocline is chosen as constant, that is 
po = 0.2 g crrT 3 from the standard solar model. The viscous heat- 
ing rate due to the MRI-driven turbulence is found to reach up to 
0(1) erg cirT 3 sec" 1 with plausible physical parameters. 

The typical traveling time T travc i of the poleward advecting 
matter due to the turbulent angular momentum transport is esti- 
mated by the diffusion time, that is T travcl =s l 2 /v turb , where / is the 
typical size of the MRI-active region measured along the latitude. 
Then the total energy gain \p of the advecting matter through the 
MRI-active region can be evaluated as, without depending on the 
arbitrary parameter /, 



</* = ^travel = PoM^ 2 / 



(6) 



where Sl tx = 2.69 x 10~ 6 sec -1 is the angular velocity at the 
tachocline. The typical size of the MRI-active region is chosen as 
'active = 2n x 0.7R Q /& here. The poleward advection flow is thus 
heated up by AT = if//(pk B /m) =s O(l0) K during which it goes 
through the MRI-active region around the tachocline. 

We would like to stress that the typical timescale of the pole- 
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Figure 3. Turbulent viscous heating rate e sustained by the MRI-driven 
turbulence as a function of co-latitude 8 in the cases with / = 0.01, 0.05, 
and 0.1. The region with the turbulent heat production closely overlaps the 
area indicating the entropy rise required from the thermal wind balance. 



ward flow traveling the MRI active region becomes 



T travcl - 3.2 x 10 



sec 

-l 



(7) 
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This is almost comparable to one solar cycle. Then we can predict 
that the poleward flow resulting from the MRI-driven turbulence 
has a propagating velocity 0(1 00) Cm SeC (— /active /"^travel) lit the 
higher-latitude tachocline. The physical picture we have arrived at 
here is more or less speculative, but provides new insight into the 
role of the MHD turbulence driven by the MRI in the solar interior. 
Further numerical study on nonlinear properties of the MRI around 
the tachocline should be prompted to make this picture robust. 

Based on the current understanding, roughly the upper third 
of the tachocline is the overshoot layer and the remaining lower 
two-thirds have the same thermal properties as the radiative zone 
(c.f., Gilman 2000; Miesch et al. 2005). It would be thus expected 
that the MRI-driven turbulence nonlinearly interact with the con- 
vective overshoot (i.e., overshooting plumes) in the upper third of 
the tachocline region. 

The impact of the convective hydrodynamic turbulence on the 
MRI-driven turbulence is studied by Workman & Armitage (2008). 
They find that the convective turbulence modifies the properties of 
the MRI-driven turbulence with the effect varying depending on the 
spatial scale of hydrodynamic forcing. Large-scale forcing surpris- 
ingly boosts the intensity of the MHD turbulence with leaving the 
character of the MRI unchanged. Low amplitude small-scale forc- 
ing may modestly suppress the MRI. 

This implies that the nonlinear properties of the MRI in the 
overshoot layer also depend on the characteristics of the overshoot- 
ing plumes (such as, overshooting velocity, its latitude-dependence 
and/or filling factor of plumes which is the volume occupancy of 
the turbulent plume in the total volume of the system) although 
they still remain unsettled in the solar interior because the numeri- 
cal simulation has not yet fully resolved the energy dissipation scale 
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(Gilman & Fox 1997; Brummell et al. 2002; Rempel 2004). The de- 
tailed numerical study of MRI-driven turbulence impacted by the 
turbulent overshooting is beyond the scope of this paper, but will 
be a target of our future work. 



4 SUMMARY 

The MRI active venue is clarified on the basis of the linear the- 
ory coupling with the observed rotation profile and the standard 
solar model. It is found that the MRI activity is confined to the 
higher-latitude tachocline and the lower-latitude near-surface shear 
layer. It is especially interesting that the MRI active tachocline re- 
gion closely overlaps with the area indicating the steep entropy rise 
which is required from the thermal wind balance in the solar inte- 
rior. 

The MRI-driven turbulence would play a crucial role in 
achieving the sun's thermal wind balance by sustaining the excep- 
tional heating and equatorward angular momentum transport. The 
warm pole with AT =s (3(10) K would be the natural outcome of 
the MRI energized turbulent activities in the tachocline region. 

The impact of the latitudinal entropy inhomogeneity on the so- 
lar rotation profile has been investigated numerically by adding the 
mechanical/thermal forcing effect on the bottom of the convective 
zone (Rempel 2005; Miesch et al. 2006). These commonly indi- 
cate that the warm pole in the tachocline is key to explaining the 
rotational features of the sun. The MRI energized turbulence could 
maintain the warm pole expected from the thermal wind balance 
and the numerical studies. 

Recently, Balbus et al. (2009) has found that, by assuming 
the coincidence between the isorotation surface and isocontour of 
6S = S - S', we can obtain the exact solution of the thermal wind 
equation (3} by the characteristic method, where S and S' are re- 
garded as the entropy distribution with and without the rotational 
effects. We stress that the characteristic curves they obtained can 
correctly reproduce the rotation profile of the mid-convective zone 
(see also Balbus 2009). This implies that, at the quasi-steady state, 
the baroclinic effect has a main role in sustaining the observed ro- 
tation profile. 

The role of the baroclinic effect is fully controversial at the 
current status. Although the importance of the Reynolds stress in 
the convective zone is indicated by large scale simulations (e.g., 
Brun et al. 2004), it does not reach consensus yet because the small- 
est size of the convective cell and the turbulence cascading process 
are inevitably restricted by the grid size even with the current huge 
computational facilities. The global MHD simulation which can 
capture the localized MRI growth and the other local scale phe- 
nomena will help us to correctly understand the physical origin of 
the solar rotation profile and the related magnetic dynamo process 
in the future. 

The rotation profile of the sun reflects much physical informa- 
tion of the solar interior. The helioseismic study covering the sun's 
polar region will reveal further details of the internal flow proper- 
ties and the roles of the MRI in the thermal wind balance of the 
solar interior. 
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